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A recently demonstrated quantum electron pump is discussed within the framework of photon- 
assisted tunneling. Due to lack of time-reversal symmetry, different results are obtained for the pump 
current depending on whether or not final-state Pauli blocking factors are used when describing the 
tunneling process. Whilst in both cases the current depends quadratically on the driving amplitude 
for moderate pumping, a marked difference is predicted for the temperature dependence. With 
blocking factors the pump current decreases roughly linearly with temperature until ksT w hui is 
reached, whereas without them it is unaffected by temperature, indicating that the entire Fermi sea 
participates in the electronic transport. 



Tunneling of electrons through classically forbidden re- 
gions is one of the major paradigms of quantum mechan- 
ics. Although our understanding has greatly improved 
over the past decades, some vital aspects of this pro- 
cess are still subject to fierce debates. For example, 
no consensus has been reached at yet whether so-called 
final-state blocking factors exist in the tunneling pro- 
cess across a barrier sandwiched between two conductors 
(Fig. Q), to enforce Pauli's exclusion principle that no 
two electrons may occupy the same quantum state. This 
seemingly innocent question is related to the question 
where in the so-called Fermi sea of conducting electrons 
the current actually flows, and thus has far-reaching con- 
sequences for our understanding of quantum transport. 
There are two schools on how to calculate the tunnel cur- 
rent, one that insists on using the blocking factors, and 
another that rejects them [Q. The dilemma is that both 
schools almost always seem to give identical answers for 
the current. In this Letter we point out that this result 
is fundamentally related to time-reversal symmetry, and 
study a generic system where this symmetry is broken to 
gain insight into the nature of blocking factors. 

Consider a tunneling barrier impeding the current flow 
as depicted in Fig. [l| According to the first school, we 




FIG. 1: Single-barrier scattering state with transmitted and 
reflected amplitudes, populated according to a distribution 
function / L characterising the left side. In the Pauli pic- 
ture, the tunneling probability incurs an additional final-state 
blocking factor 1 — / R on the right. 



start by calculating the coherent scattering states of the 
system consisting of an incident plane wave, a reflected 
wave, and a transmitted wave on the far side. These 
scattering states travel either from left to right or the 
opposite way. In the stationary limit they can simply be 
populated according to distribution functions / L and / R 
governing the asymptotic contact regions on either side. 
This yields for the current jl) 

Is = ^ J dEdE'V ± T+(E',E)f L (E) 

- ^ J olEdE 1 V 1 _T-{E,E')f K {E'), (1) 

where T>± is a density-of-states factor T + (E',E) is 
the transmission probability for scattering states incident 
from the left at energy E and emerging at the right at E' 
(7^ E in general), and T~ is defined in a similar manner 
for the reverse direction. According to Eq. ([|) all scat- 
tering states contribute to the current, even those far 
below the Fermi surface — as recently discussed for the 
Fermi pump |3] . It is only when the transmission proba- 
bilities are symmetric, T + (E',E) = T~(E',E), that the 
net current seems to stem from electrons in the immedi- 
ate vicinity of the Fermi surface only, as in this case all 
electrons further down cancel each other. 

An alternative and widely used recipe for calculating 
the tunneling current in Fig. ^ is based on the transfer- 
Hamiltonian formalism, originally put forward by Op- 
penheimer and later refined by Bardeen and others 
Q. Here the system is split into two subsystems, the left- 
and the right-hand side, with a common overlap in the 
central barrier region, and the result for the current can 
be expressed as 

2 ' f dEdE' T> ± T+(E', E)\l - f R {E')]f L (E) 



h 
2e 

77 



dEdE' X> ± T- (£,£')[! - f L (E)]f R {E'l2) 



* wagner@phy.cam.ac.uk 



which differs from Eq. (Q) by the Pauli blocking factors 
1 — /. These factors are introduced — more or less ad 
hoc — using the intuitive argument that an electron tun- 
neling from one side to the other needs to find an empty 
final state on the far side to tunnel into [pf. 
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FIG. 2: Quantum-dot pump with broken time-reversal sym- 
metry. [A] Schematic drawing; [B] simple ID model based on 
driven 5-function barriers. 



Clearly, the philosophies underlying Eqs. (|l|) and (g) 
are entirely different: For the former we assume that the 
left- and right-hand side of the barrier belong to the same 
system and that the scattering state extends coherently 
across the barrier, whereas in the latter we consider real 
incoherent transitions between initial and final states on 
opposite sides of the barrier belonging to different sub- 
systems Q. Both views have their merits. Certainly in 
the limit of a very transparent or even vanishing barrier, 
one has to regard both sides of the barrier as part of the 
same system, in which case there is no place for blocking 
factors. But equally well, if tunneling is weak, one can ar- 
gue that the two sides do form two different subsystems, 
and that final-state blocking factors are mandatory to 
guarantee Pauli's exclusion principle and to prevent an 
"overflow" of the final state. 

Crucially, this analysis also applies to a system with lo- 
cal time-dependent driving forces as its leads, and hence 
the distribution of incident electrons, are not affected by 
the driving field. For a system with time-reversal symme- 
try, it can be shown that the transmission probabilities 
are such that any channel E <-> E' is traversed with the 
same probability in either direction, i.e. that microre- 
versibility holds, T+(E',E) = T-{E,E'). Under this 
condition the cross terms / L / R arising from the Pauli 
blocking factors cancel and Eqs. (|l|) and (||) yield iden- 
tical answers for the total current, which is the quantity 
usually measured in experiment. It is probably because 
of this indecisive result that both schools have coexisted 
for so long. The only difference between the two is their 
prediction where in the Fermi sea the current flows: With 
blocking factors, the current is forced to flow close to the 
Fermi surface, as in this picture lower-lying final states 
are blocked, whilst in the formulation based on scattering 
states the current flows in the entire Fermi sea. 

However, if time-reversal symmetry does not hold, 
Eqs. (0) and (|) yield different answers even for the mag- 
nitude of the current itself as the cross terms do not can- 
cel any longer. This is very intriguing as it gives us hope 
to come closer to an experimentally verifiable answer re- 
garding the existence of Pauli blocking factors. 

Our analysis is inspired by a recent experiment by the 
Marcus group B. A semiconductor quantum dot with 
source and drain point contacts has two additional lat- 
eral gates to which ac voltages of relative phase <f> are 



applied (see Fig. [2] A). Time-reversal symmetry is broken 
unless <fi is an integer multiple of it. Averaging over differ- 
ent dot configurations, they measured fluctuations of the 
emf voltage generated between the source and drain con- 
tacts. Previous theoretical studies employed the concept 
of adiabatic pumping |j) . We take a complementary per- 
spective and view the pump current (which exists even 
in the absence of any applied source-drain bias) as due 
to photon-assisted tunneling Q . 

As a model system of a driven dot strongly coupled 
to its leads we consider two harmonically oscillating 5- 
function barriers of equal strength U ac a distance d apart 
in a ID potential, with a variable phase difference <fi in 
the ac signals as depicted in Fig. ^ B ]rj{ |. This is a 
very simplified model of the Marcus pump, but never- 
theless it turns out to exhibit many of its characteristic 
properties. For calculating the transmission probabili- 
ties T(E',E) across the dot we take advantage of the 
fact that we can split this problem into two parts: If we 
know the transmission and reflection amplitudes for each 
barrier separately, we can use the Fabry-Perot method of 
raytracing known in optics to calculate the interference 
pattern due to multiple reflections between the barriers. 
As in optics the partial interference amplitudes can be 
summed up to all orders in a geometric series, yielding 
for the transmission amplitudes at the far side |l2| 
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where K = Qr L Qr R describes one full round trip of the 
electron between the two barriers L and R, starting and 
ending at the R barrier. Each round trip an electron 
at energy E n = E + nfru picks up a phase factor Q n — 
exp[i(k n d+6)] consisting of two parts: The phase k n d in- 
curred after travelling a distance d with wave vector k n = 
\/2mE n /Ti, and a fixed but unknown phase 6. This latter 
phase is introduced to account for the random changes 
in magnetic field and dot geometry employed in the ex- 
periment to perform ensemble averages. In general, the 
phase 9 depends on the electron energy, but for simplic- 
ity we ignore this and assume 8 to be equally distributed. 
Under this condition ensemble averaging simply means 
averaging over 9 at the end of the calculation. 

Due to the discreteness of the photon energy, electrons 
emerge on the far side of the barriers at energies differing 
from their original energy E by multiples of Tim: E' = E 
+ nhw, the so-called sidebands. Consequently, all ampli- 
tudes in Eq. (^|) have to be interpreted as matrices, and 
the transmission probability takes the form T(E', E) — 
f n {E)S(E + nhw-E'). 

Before presenting numerical results for strong driving, 
it is instructive to study the weak-driving limit, which 
can be solved analytically. Defining T c (k, if) = ij 2 (fco/4fe) 
cos 2 [d(ko + k)/2 + if/2] where i] = V ac k /E is the di- 
mensionless pump amplitude and kg the wave vector 
of the incident electron, we obtain for the transmission 
probabilities in the first three sidebands up to rj 2 |13 
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FIG. 3: Dependence of (if) 1 / 2 on the phase shift <f> for weak 
and strong driving fields. Parameters: d — 0.2 /im, f = 
10 MHz, m = 0.067 m , E F = 12meV, T = 0.1K. 
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FIG. 4: (7j(7r/2)) 1/2 and {Vi{n/2)) 1/2 , both based on Eq. (@), 
as a function of the driving strength 14c for a range of tem- 
peratures. 



-Te(fci,20-0)-T e (fc_i,20 + 0) 
£+(0,0) = -T c (-fc ±1 , T 0). (4) 

The corresponding transmission probabilities T~ for the 
reverse direction are obtained by substituting (0, 0) — > 
(—</>,#). As expected, the net transmission probability 
T no t = T + - T~ vanishes for = ri7r when time-reversal 
symmetry holds, and is maximal at 4> = tt/2. In agree- 
ment with experimental findings, the pump current at 
zero dc bias, being proportional to T nct , scales with rj 2 
for weak driving. In the experiment the photon energy 
is 6 orders of magnitude smaller than the Fermi energy. 
Expanding Eq. (H) in this limit the current turns out to 
be linear in the driving frequency, in perfect agreement 
with the experimental results available | Q . 

After integrating over 9 the ensemble average (T net ) 
is found to be identically zero up to order rj 2 . This im- 
plies that the direct current (I) is orders of magnitude 
smaller than its fluctuations for all but the strongest driv- 
ing fields, and hence we will concentrate on the mean 
square average (J 2 ) 1 / 2 to study the fluctuations instead 
p5[ . Figure || shows the dependence of (if) 1 ^ 2 on the 
phase shift </>, where we have used Eq. (|^) with an adap- 
tive number of photon sidebands to determine the full 
non-linear transmission probability, and Eq. (0) to finally 
calculate the current. For small driving amplitudes we 
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FIG. 5: (^(7r/2)) 1/2 , based on Eq. (|), as a function of the 
driving strength 14c for a range of temperatures. 



find a sin0 behaviour in agreement with Eq. (Q), whilst 
for stronger driving nonharmonic features appear, both 
of which are in qualitative agreement with experiment. 
A similar behaviour is observed when using the formula 
(||) based on Pauli blocking factors instead, except that 
in this case the pump current is much smaller. 

For the remainder of this Letter we will fix the ac phase 
shift (f) to 7r/2, the point of maximal time-reversal asym- 
metry. The current fluctuations (/| (tt/2)) 1 / 2 calculated 
using the formula (|l|) without blocking factors are illus- 
trated in Fig. [| as a function of i], taken at the Fermi 
surface. For small driving amplitudes up to rj « 0.05 the 
fluctuations increase quadratically with r\ as suggested 
by Eq. before eventually starting to decrease. In 
the Marcus experiment, rather than measuring the pump 
current, the emf voltage generated was studied. Defin- 
ing the emf voltage as the difference in chemical poten- 
tials between the left and right leads necessary to make 
the pump current vanish, we find that its fluctuations, 
also shown in Fig. ^ (right-hand axis), have virtually the 
same dependence on r\ as the current fluctuations. Both 
the quadratic rise as well as the levelling off for stronger 
pumping have been observed in experiment. The maxi- 
mal emf fluctuations generated, 15 nV, is only one or two 
orders smaller than in the experiment, which is a quite 
reasonable agreement given the simplicity of the model. 

The corresponding results based on the alternative for- 
mula (||) for the current which does include the Pauli 
blocking factors are illustrated in Fig. g. Similar to the 
case without blocking factors of Fig. [l| the current fluctu- 
ations (Ip(Tr/2)) 1 / 2 (as well as the emf fluctuations, not 
shown) display an rj 2 behaviour first before eventually 
starting to saturate. Yet, not only do the magnitudes of 
these results differ substantially from Fig. [I| for the same 
driving amplitude, but now there is a very pronounced 
temperature dependence as well! 

This temperature dependence is worked out in more 
detail in Fig. ^ using a fixed driving strength of rj = 0.47. 
Whilst without Pauli blocking factors the variation with 
temperature is minimal, an almost linear dependence is 
observed when they are included. This distinctly dif- 
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ferent behaviour in the low-temperature regime is most 
easily understood when looking at the number of states 
in phase space effectively available for transport: The 
blocking factors force the current to flow within a few 
k^T of the Fermi surface (see Fig. |l|), and as the tem- 
perature approaches zero, this range of active current- 
carrying states eventually diminishes to a minimal width 
of a few fouj, which in the experiment is much smaller 
than k^T . However, since each state can only carry a 
certain maximal load, it follows that the pump current 
must also decrease with temperature, until it settles for 
a residual value once k^T « fiu is reached — with our 
parameters at ~ 0.5 mK. On the other hand, without 
blocking factors there are no phase-space restrictions, in 
which case the pump current flows in the entire Fermi 
sea and thus is largely immune to changes at the Fermi 
surface brought about by temperature. 

In the Marcus experiment the pump current is found 
to increase when lowering the temperature, and appears 
to level off at w 0.1K, where phase- breaking events be- 
come less important H. Such events are not included 
in our theory, and we can therefore not expect to re- 
produce the high-temperature behaviour. However, the 
experimental finding of a saturated pump current in the 
low-temperature limit is (if genuine and not due to ther- 
mal decoupling) only consistent with our results based on 
the scattering-state approach, Eq. (fy), but not with the 



formulation relying on Pauli blocking factors, i.e. Eq. (||). 

Being able to prove the (non)existence of Pauli block- 
ing factors has drastic consequences for deciding whether 
the current flows in the entire Fermi sea or at its surface 
only. In this Letter we have demonstrated that a pow- 
erful tool for studying this issue is to look at the tem- 
perature dependence of the pump current when breaking 
time-reversal symmetry. Although our model system is 
simple, the conclusions about the low-temperature be- 
haviour, being drawn from phase-space considerations, 
are clearly of a much more universal nature. 

The author acknowledges discussions with C. Marcus, 
and support by the EU (FMRX-CT98-0180). 
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FIG. 6: Temperature dependence of (I 2 (n /2)} 1//2 with and 
without Pauli blocking factors. 
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